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Note on the Theory of Simultaneous Linear Differential 
or Difference Equations with Constant Coefficients. 

Bf J. J. Sylvester. 



This theory is virtually the same for differential as for finite-difference 
equations. The mere verbal part of the exposition being somewhat easier for 
the former of the two, I shall prefer in the first instance to deal with them, 
although the applications are more interesting when made to bear on the latter. 
Simple to the last degree as are the method of solution and the nature of the 
result, I do not find the one or the other set out, or even indicated, except in "the 
most perfunctory manner, in the ordinary text-books. This brief notice, designed 
for the junior readers of the Journal, is intended to supply the lacuna. 

Let u^ k denote a linear function, with constant coefficients, of u k and of its 
first Sj derivatives in respect to L 

Let «i, i + «i,»+ • • • ■ +«i,i = 

«s,i + **»,» + • • • • +%,t=0 



be the system of differential equations proposed for integration. 

Call ex + e 2 -f- . . . . + e t = a . 

The process of arriving at the reducing equation for any one of the variables 
is after the manner of the dialytic method of elimination, viz.: 

Along with the first equation take each of its (a — ex) tb derivatives, with 
the second equation each of its {a — e 2 ) th derivatives, . . . and with the * th equa- 
tion each of its (cr — e^) th derivatives. 

There will thus come into existence (* + 1) a — fa + e 2 . . . . + £«) *• e. 
i (c + 1) — cr equations between the i (<r -f- 1), quantities 

Ox, 8 t G>x, .... 5f% 

S t o 2 , . . . ._ 8fa 



CO, 



COi, § t Q t , .... StO t 
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If we omit those which appear in any one of the lines above written, there 
will remain (a + 1) (i — 1) or i (a + 1) — a — 1 which might be eliminated 
between the i (<r + 1) — a equations, and there would thus result an equation 
between the quantities contained in the omitted line. This elimination, it 
will presently be seen, there is no occasion to perform ; the noticeable algebraical 
fact about it is, that supposing it were performed, the form of the equation 
resulting between a k , 8 t a k , . . . $a k is invariable, whichever of the numbers 
1, 2, 3, . . . i be the value assigned to h. 

Let the order of the highest derivative of each a be reduced by one unit 
below the highest order previously taken, then there will be ia — cor (i — 1) a 
equations connecting the ia quantities 

a> 1} $ t o lt .... $"- 1 o 1 
Oi, £ t a 2 , . . . . ^ -1 % 



>a lt . . . . tj ] g> 4 



and accordingly, if we omit the a quantities which appear in any one (say the 
first) of the above lines, the remaining (i — 1) a quantities may each of them be 
expressed as linear functions of % and its (<r — 1) derivatives : but the elimina- 
tion previously indicated would lead to a homogeneous linear equation between 
% and its a derivatives, and if in that, each argument hfOi be replaced by h K 
and \, \, . . . % a be the a roots of the algebraical equation so formed, it follows 
from the ordinary theory for a single equation that ^ (provided the given equa- 
tions, and consequently the resulting ones, be left in their general form) will be of 
the form 1 C 1 e h S + 1 G 9 e^* t + . . . -f- 1 C a .e h ° t , and consequently by virtue of the pre- 
vious remark o 3 , o 3 , . . . a Jc will be of the same form as % (but, of course, with 
different coefficients), that is to say, the a roots h lf \, . . . h a are the same for 
the equation in a k as for the equation in a 1 , so that the coefficients in the equa- 
tion between % and its a derivatives are, as premised, independent of the value 
of It,. 

Finally, to determine the equation whose roots are \, h 2 , . . . h„, let 
1 Ce M , one of the terms in the general value, be taken as a particular value of %, 
which with corresponding values of the other o's will serve to satisfy the given 
equations ; %, a 3 , . . . a t being each of them linear functions of a and deriva- 
tives of a, must be of the forms *Ce M , s Ce ht , . . . W, so that o 1( %, . . . o< 
and the derivatives of each of them will contain the -common factor e' rf , and by 
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substitution in the original equations we shall obtain a system of simultaneous 
algebraical equations leading to the equation 

-"l, 1) -"1, 2 • • • • Ri, s 

R», 1> R%, 2 • • • • -*H, i 







Ri, l! ^i, 2 • • • • Ri, i 

where in general R pq is what u Pt q becomes on writing A" in place of %u q . 

The above determinant of the i th order will be of degree e x + e 2 . . . + e* , 
*. e. of the degree <r (for the general case) in h , and the roots of the equation 
will give the a values 7^ 1 , h 2 , . . . h a . 

It follows, therefore, that the result of the hypothetical elimination in the 
first instance referred to will be a linear function of 5*%, h" t ~ 1 %, . . . S t a k , o k 
of which the Coefficients will be identical with the coefficients of h", ~h°~ x ,-. . . h, 1 
in the above determinant. Hence no matter now what special values may be 
attributed to the coefficients of the given equations, the result last obtained 
•remains of universal validity — without excepting those cases in which the result 
of the hypothetical elimination would be such that the corresponding algebraical 
equation possess equal roots, although in those cases the form assumed in the 
course of the argument for the value of a x (viz. a linear function of exponen- 
tials) ceases to hold good. Neither for the same reason need any exception be 
made for those cases where the number of terms in the equation to o k falls 
below a on account of one or more of the leading coefficients in the result of 
the hypothetical elimination becoming zero : the degree to which h rises in the 
determinant will be in all cases the right degree, whether it reaches the extreme 
possible limit a or falls below it. 

The result obtained may be briefly summarized as follows. 

if (M) » + (M)y + + (M) » = ° 

0M.)« + (M)y + + CM) * = o 



(%$*) x + {aA)y + + (»A) 2 = 

(each <?> , 4 1 > • • • j a standing for a rational- integral functional form) then will 

(R8 t ) x=0 (RS t ) y = . . . . (Rh t ) e = 0, 
where R(8 t ) is the resultant in respect to x, y, . . . z of what the above 
equations become when S t is treated as an ordinary algebraical quantity ; under 
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which form the proposition (by virtue of Euler's method of multipliers) becomes 
so nearly intuitive as to abrogate all necessity for any other demonstration.* 

To pass to the parallel and more important theory in finite differences, it is 
only necessary to interpret u Jt ft to signify a linear function, with constant coeffi- 
cients, of (a k ) t , (%)(4-i, . . . (o A ) <+ ,.,, where t is the integer independent variable, 
(say (a k ) t and its e, difference-augmentatives), and instead of taking the differ- 
ential derivatives of any one of the given equations, to take the corresponding 
difference-augmentatives. Then by precisely the same reasoning as before we 
shall have 

(>it+*-\- Ba t+a _ 1 + . . . .-\-La t =0, 
B, G, . . . L being so taken as that h° + Bh°~ x + . . . + L shall be the 
determinant represented by the same form of matrix expressed by B's as before, 
but where JKp i9 is obtained from u p ^ q by writing h 9 in lieu of any argument 
&t + Q which occurs in it. 

The simplest example that can be given is where i = 2 , e 2 = e 2 = 1 , 
«i, i= — m+i + a *lt> **i, a=&0«> 



u. 



'2, 



x — GYlt Ui, %=— Ot + l + dd t i 







this was the case which occurred in the article on the extension of Tchebycheff 's 
theorem, in the last number of the Journal, leading to the equation 

a — /I b 
c d — % 
and to expressions for v\ t , t as linear functions of h\ , X\. 

It may also be remarked that this same case gives an instantaneous solution 
of the problem proposed and successfully treated by Babbage in his Calculus of 
Functions, more than half a century ago, and since revived in connection with the 
theory of substitutions (Serret, Alg. JSup. 4 ed. torn. 2, pp. 256-262). The problem 

is to find tyx = . , so that <p*x , say -^ — ' , , shall equal x for a given value of i . 

px -J- O pi X -f- 0{ 

To find in general fix it is only necessary to solve the difference equations 

u t = au i _ 1 + «.« i _ 1 

and then ii t , v t will, if u = 1 , v = 0, coincide with a t , (3 it and if u = , v = 1 
with a t , \ . 

* I regret that this simple reflection did hot present itself to my mind before the preceding investi- 
gation, the necessity for which it does away with, had been set up in print. It of course applies equally 
well to the analogous proposition for finite-difference equations (w s , v t , . . . . being substituted for 
x , y . . . . , and 1 + A for 6,) . This last named proposition, limited to the case of equations of the first 
order, is the foundation-stone of my new theory of Matrices regarded as Quantities, i. e. as subject to 
every kind of functional operation which ordinary arithmetical or algebraical quantities are or can be 
subject to : but though so important and so easily established, I know not where it can be found expli- 
citly stated. 
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= 0, 



Thus calling o x , o % the two roots of 

p + a a 

P —p + b 

a t will be of the form C (p{ — pi) and {3 t of the same form except as to C , say 
T (p{ — p|) . Also a { , bi will be of the forms Qpl + <7 2 p| , I\pi + T%pl , where 
C\ + @% — 1 > r x + r 2 = 1 , and the required condition will be fulfilled, pro- 
vided only that pj = p| , or say 

.„ / Xn , , . )jc \ 

p 1 = K I cos — — hv — 1 sin — — ) 

t-t / to . . Xn \ 

p 2 = K I cos — V — 1 sin -r- ) 

*". e. if (a + by — 4 (a& — a/3) f cos — — J = 1 , X having any integer value (which 

without loss of generality may be taken inferior to i) except zero.* 

If X = , the two roots of the equation in p become equal and the form of 
the solution changes into 

When u = 1 and « = then % = a , v = /3 , 

Ci=l, Oi' = 0, '(7,= -f-l, Ci'=4. 

r I 

and when rf = , v = 1 , % = a , v g = b , 

C, = j, C' = j~l, C 2 = 0, C 2 '=l, 

and <p l x = ■ .7" , ~r — s • which cannot be periodic for any value of * , 

r ifix -+- p -j- (6 — ,o) * r ^ 

,, . , (a — p) x 4~ a a — p a . a — b 2a 

and when % = co becomes -^ — —^ — ' — = - = - , *. e. = „. or =- , 

fix + 6 — ,o /* 6 — p ' 2,3 a — 6 

so that c^tc in this case continually converges to a constant limit. 

I may add that tfx converges to a constant limit not merely when the roots 

Pi i pa of 

a — p a 

P &-p 

are equal, but whenever they are real. For the general form of q> l x , it may 

easily be found, is 

\_{p% — a) Pi + (pi — a) gl] x + a (p{ — p\) 
g {pi — pi) x + [(ft. — b)pj + (Pi — b) pj] 

* There will thus be (i — 1) values of % which will each give a distinct admissible solution of the 
problem of periodicity, but of course only those values of A which are relatively prime to i will give 
primitive solutions. If i — i'$ the effect of making X^zVS will be to make fx — x by virtue of itsr 
making f'x — . 
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which if p 2 *> pi when i = ao becomes -^_ — I =- = — —^ or -. where 

1 * ^ > 1 — p®-\-pi — j3 b — p x 

p x signifies the smaller of the two roots o x , p 2 5 or m other words when 

* r (IT I (7 

a — 5>2V«|3, the limiting value to <p*x, when $x represents J~ : > is 

v a— j+v(«— ) «p ^ with the understanding that the quantity under the radical 

2/5 

sign is to be taken positive. 

So, if 

x i+1 : ij i+1 : z i+1 = axi + by t + <w, : a'asj + %* + c ' 2 j : a '% + 6 "2/* + c " z i > 
when all the roots of the determinant 

a — X b g 

a' V— X d 

a!' b" c"— X 

are real, the point x i} y t , z t , as i increases, will be found to approach indefinitely- 
near to a fixed straight line ; and if all the roots are equal, to a fixed point. 

The condition of the system of ratios x t : y t : z t being periodic and having a 
period m is tantamount to the condition that the ?n th power of the matrix 



shall be the matrix 



a 


b 


e 


a' 


V 


d 


a!' 


b" 


d< 


1 











1 











1. 



The complete solution of this problem, and of the more general one of ex- 
tracting the wi th root of any unit-matrix (*'. e. a matrix in which each element in 
the principal diagonal is unity, and the rest zero), which constitutes the ultimate 
generalization of Babbage's problem and is soluble by the same method, will 
probably appear in a memoir on matrices, in the forthcoming number of the 
Journal. 

In general, for a matrix of the order a , the number of m th roots is m!* and 
each of them is perfectly determinate. But when the matrix is a unit-matrix 
or a zero-matrix (the latter meaning one in which every element is zero) 
there are distinct genera and species of such roots, and every species contains 
its own appropriate number of arbitrary constants. 



